Pairing in the framework of the unitary correlation operator method (UCOM): 

Hartree-Fock-Bogoliubov calculations 
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In this first in a series of articles, we apply effective interactions derived by the Unitary Correlation 
Operator Method (UCOM) to the description of open-shell nuclei, using a self-consistent Hartree- 
Fock-Bogoliubov framework to account for pairing correlations. To disentangle the particle-hole and 
particle-particle channels and assess the pairing properties of Vucom, we consider hybrid calculations 
using the phenomenological Gogny D1S interaction to derive the particle-hole mean field. In the 
main part of this article, we perform calculations of the tin isotopic chain using Vucom in both the 
particle-hole and particle-particle channels. We study the interplay of both channels, and discuss 
the impact of non-central and non-local terms in realistic interactions as well as the frequently used 
restriction of pairing interactions to the 1 5o partial wave. The treatment of the center-of-mass 
motion and its effect on theoretical pairing gaps is assessed independently of the used interactions. 

PACS numbers: 21.60.-n,21.30.Fe,13.75.Cs 



I. INTRODUCTION 

Recent years have seen a revival of nuclear structure 
physics, motivated by new experimental advances in the 
use of radioactive beams in existing and proposed facil- 
ities at GSI/FAIR, RIKEN, GANIL, and other labora- 
tories worldwide, as well as new theoretical approaches 
to the nuclear many-body problem. The application of 
effective field theory (EFT) and renormalization group 
(RG) methods has provided systematic approaches to the 
construction of effective nuclear interactions that main- 
tain a stringent link to QCD, either directly on the formal 
level or by the reproduction of low-energy observables like 
NN phase shifts and deuteron properties. The former 
have culminated in the derivation of a consistent set of 
two- and higher many-nucleon interactions in the frame- 
work of chiral EFT at next-to-next-to-next-to-lcading or- 
der (N3LO) while the latter have revealed the uni- 
versal aspects of realistic NN interactions by decoupling 
low- and high-momentum modes via RG decimations in 
the case of Vjow-k > or unitary transformations in the 
Similarity Renormalization Group (SRG) Q. 

While starting from a different premise, i.e., the ex- 
plicit treatment of correlations induced by the repulsive 
core and the tensor force of realistic NN interactions like 
Argonne V18 or CD-Bonn (see @ for a review) , the Uni- 
tary Correlation Operator Method (UCOM) [1,0] shares 
many characteristics of the RG-derived low-momentum 
interactions. This is particularly true for the SRG ap- 
proach, where the dynamical generator of the unitary 
transformation is related to the generators of the UCOM 
transformation [1,0]. 

In our previous works, correlated interactions derived 
in the UCOM framework, referred to as Vucom in the 
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following, have proven their merit in a wide range of 
many-body methods, from ab-initio calculations of light 
nuclei in the No Core Shell Model 0,1111, HI to Hartree- 
Fock (HF) and HF-based extensions like Many-Body Per- 
turbation Theory 0, H], RPA 0, and Second 
RPA fir) . Since HF-based approaches do not account 
for pairing correlations, they are expected to work best 
for closed-shell nuclei. The aim of this first in a series 
of papers is the extension of our calculations to open- 
shell nuclei by constructing a self-consistent Hartree- 
Fock-Bogoliubov (HFB) framework (see, e.g., (l8j). 

Initial attempts to solve the HFB equations using (at 
the time) "realistic" interactions [HI, [2(j were hampered 
by the strong repulsive core of the NN interaction, which 
leads to infinities in mean-field methods because Slater 
determinants arc inherently unable to describe the re- 
quired correlations. Brueckncr's G-matrix approach plj ] 
provided a way to deal with this problem by resum- 
ming particle-particle ladder diagrams, leading to a well- 
behaved "tamed" interaction, but remained problematic, 
e.g., due to the starting-energy dependence. At the same 
time, it was observed that a similar resummation was 
not required to regularize the NN gap equation (see, e.g., 
[221 , [23j). and, in modern terms, the "bare" interaction 
could be used directly, prompting a disparate treatment 
of the particle-hole and particle-particle channels in self- 
consistent field calculations in the following decades. 

Parallel to these original Hamiltonian-based ap- 
proaches, Negele and Vautherin introduced the Density 
Matrix Expansion 24[ in nuclear physics. While their 
work was tied to the Hamiltonian-based approaches by 
using similar concepts as in G-matrix methods, it also 
provided a foundation for the form of Skyrme-type en- 
ergy functionals, and paved the way for phenomenolog- 
ical Density Functional Theory (DFT), which became 
the standard framework for self-consistent field methods 
until today (see, e.g., [25| for a comprehensive review). 
While current phenomenological density functionals are 
able to describe nuclear bulk properties like binding en- 
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crgics and charge radii with high accuracy near the val- 
ley of stability, they often perform inadequately in the 
description of spectroscopic observables or exotic nuclei. 
As a result, considerable theoretical effort is under way 
to improve the phenomenological functionals (see, e.g., 
[26l . l27l|). or to construct it from first principles by ap- 
plying EFT methods (28J. Given the guiding principles 
of EFTs with respect to consistency, one then has to ask 
whether one should actually demand the treatment of the 
particle-hole and particle-particle channels on the same 
footing — especially since effective interactions derived 
in the (S)RG or the UCOM approaches no longer require 
resummations in the particle-hole channel. 

A strong argument in favor of such consistency was 
encountered in phenomenological DFT in recent years. 
Many-body methods beyond the mean field, e.g., Gener- 
ator Coordinate approaches, involve configuration mixing 
of non-orthogonal Slater determinants, and rely on sub- 
tle cancellations between singular terms in the particle- 
hole and particle-particle channel (see, e.g., [29l. |30|]). A 
Hamiltonian provides an ideal starting point for these 
methods, because the use of the same interaction in both 
channels automatically guarantees these cancellations, 
whereas one has to go to some lengths to implement a 
"regularization" scheme in DFT to remove spurious con- 
tributions to the energy and other expectation values due 
to the use of separate particle-hole and pairing function- 
als mm- 

For the reasons discussed above, the aim of this paper 
is the formulation of a fully self-consistent HFB scheme 
based on an intrinsic Hamiltonian, using a family of cor- 
related interactions based on the realistic Argonne V18 
interaction [33[ in both interaction channels. After briefly 
reviewing the basics of the HFB approach and the Uni- 
tary Correlation Operator Method in Sect. UH we dis- 
cuss certain details of our implementation in Sect. IIIII 
This includes a comparison of the convergence behavior 
of Vucom with the phenomenological Gogny D1S interac- 
tion [34| , as well as a detailed discussion of the center-of- 
mass treatment. In Sect. IIV1 we investigate the proper- 
ties of Vucom as a pairing interaction, using the Gogny 
force to generate the mean field and thereby disentan- 
gle the particle-hole and particle-particle channels. This 
section also includes a comparison with the pairing prop- 
erties of SRG-evolved interactions. Section [V] presents 
results from fully self-consistent HFB calculations with 

Vucom- 
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where Ck and c\ are annihilation and creation operators in 
the standard particle space. The Bogoliubov conditions 

U^U + V f V =1, UU t + V*V T = l, (3a) 
U T V + V T U = , UV^ + V*U T = , (3b) 

ensure that the {f3k, /?fc}fceN satisfy the canonical anti- 
commutation relations. In HFB approximation, the nu- 
clear ground state \*&) is defined (up to a unitary trans- 
formation) by the quasiparticle vacuum, 



=0. 



(4) 



Taking account of the center-of-mass contribution to 
the kinetic energy, we introduce the intrinsic kinetic en- 
ergy operator 



A 2 

~ A 2fi ' M ~ 2 ' 



where 



(5) 



(6) 



is the relative momentum and p the reduced mass, given 
in terms of the nucleon mass tun- Using T; nt , the intrinsic 
many-body Hamiltonian reads 



Hint — Tint + V — T — T c , 



V. 



(7) 



The intrinsic energy of the HFB ground state can then 
be expressed in terms of the density matrix 

Pkk ' H*|4, Cfc |*H (u*u T ) fcfc , (8) 

and the pairing tensor 



n kk , = <*| Cfe ,c fc |*) = {V*U T ) kk , 



(9) 



E[p, K ,K*] = itr(M-^r(AO, (10) 



II. FORMALISM 



A. Hartree-Fock-Bogoliubov Theory 



where the hcrmitian particle-hole field 



dE 



dpk'k \A 

qq> 



v l p qq < (11) 

kq' k' q 



The HFB approach [18( aims for a simultaneous mean- 
field description of the particle-hole and particle-particle 
channels of the NN interaction by introducing quasipar- 
ticle operators /3^}fcgN via the Bogoliubov transfor- 



and the antisymmetric pairing field 
dE 



A kk , = 



8k 



kk> 



^ 1l' J kk'qqi 



iqq' 



(12) 
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have been introduced, i and v denote the antisym- 
metrized matrix elements of the intrinsic kinetic energy 
and the NN interaction, respectively. 

The HFB ground state is obtained by performing a 
variation of the energy with respect to p and n, subject 
to the constraint 



tip = N, 



(13) 



which ensures the conservation of the mean particle num- 
ber. Carrying out the variation, one obtains the Hartree- 
Fock-Bogoliubov equations 



H 



h — X A 

-A* -h* + X 



(14) 

where we have introduced the HFB Hamiltonian TL. Due 
to the use of an intrinsic Hamiltonian, the Lagrange mul- 
tiplier A in Eq. (Tj"4")) can no longer be identified directly 
with the Fermi energy of the system (see Sect. IIIIB[) . 
Equation (fT4|) constitutes an eigenvalue problem that has 
to be solved self-consistently due to the dependence of TL 
on p and k. 



B. Spherical Symmetry 

Assuming spherical symmetry, the Bogoliubov trans- 
formation reduces to the form 

Pnljm = E ^n'n C \i'ljm + ( — ^-V^^n'n C n'lj-m i (15) 
n' 

n' 

where n is a radial quantum number and the upper in- 
dices (Ij) mark the (diagonal) angular-momentum quan- 
tum numbers. Aside from the explicit phase factor in Eq. 
(|15p . the transformation is independent of the angular- 
momentum projection m. Using the matrices U and V 
from Eq. ifPo]). one can define the reduced matrices p^ n , 
and k%$, 

Pnljm,n' I' j' rn' [VV j . 6 jj' 5n' $rnm f 

= & jf S lll S mm ^ (17) 

K n ljm,n'l'j'm' = (~ 1)^ ™ \VU T '\ J^, 5jf 5u> 5 m - m > 



(18) 



which are both symmetric and real. The antisymmetry of 
the pairing tensor is now contained entirely in the phase 
factor. Correspondingly, the reduced fields are defined 
by 



r = a a a 

-L nljm^n'l'j'm' — "jj'Ull'Uinm'*- nn ' 



(19) 



and 



^■nljm.n'l'j'm' = Sjj'Su>S m ^ m '(~l) 3 ""A®! . (20) 



With these definitions, the reduced HFB equations 
read 



hW) - A -A*') 





(21) 



C. Canonical Basis 

The canonical basis is a convenient tool for the dis- 
cussion of the HFB problem, because the HFB equations 
(TT4"]) assume the same form as in the Bardeen-Cooper- 
Schrieffer (BCS) case in this representation. It is ob- 
tained by diagonalizing the one-body density matrix l[8]). 
whose eigenvalues v 2 are interpreted as occupation prob- 
abilities of the corresponding canonical states The 
accompanying coefficients u 2 ^ are defined up to a phase 
by the condition [cf. Eq. ([3])] 



1 



(22) 



Analogous to the BCS case one can then define gener- 
alized sin gle-p article energies and state-dependent gaps 
(see, e.g., [35|) via the matrix elements 



A^ = A„„ , 
as well as the canonical quasiparticle energy 



A) 2 



(23) 
(24) 



(25) 



Here, \p) is the canonical conjugate state of (e.g., the 
time-reversed state in systems with time-reversal symme- 
try), is just the diagonal matrix element of the HFB 
Hamiltonian in the canonical basis, and therefore gen- 
erally not identical to any of the quasiparticle energies 
obtained by diagonalizing TL. In terms of the newly de- 
fined quantities, the canonical coefficients can be written 
as 



(26) 
(27) 




where we have adopted the phase conventions of [35: 



D. Gap Definitions 

Experimentally, the odd-even staggering of nuclear 
binding energies provides a clear signal of pairing correla- 
tions in the finite nucleus. This staggering is analyzed via 
the ground-state energy differences of several neighbor- 
ing nuclei. Recently, it was argued that the odd-centered 
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three-point formula 



generators 



A^(iV) 



(E(N + 1) - 2E(N) + E(N - 1)) 



(28) 

provides the clearest measure of pairing correlations 
along isotopic (or isotonic) chains, because it is least af- 
fected by particle- hole effects (see, e.g., [36|, which also 
discusses further refinements). The best way to com- 
pare these experimental "gaps" to theory would be the 
application of Eq. (|28[) to theoretical ground-state ener- 
gies. Since the treatment of odd nuclei in a HFB frame- 
work requires further approximations with respect to the 
blocking of levels by the unpaired nucleon, we defer such 
calculations to the future. We point out, however, that 
such calculations are facilitated in a Hamiltonian-based 
approach because the interaction is already completely 
determined. 

While the pairing energy in Eq. ([11])) provides an obvi- 
ous measure of pairing correlations in a theoretical calcu- 
lation, it cannot be related directly to the experimental 
gap ([28]) . Thus, one usually turns to the canonical basis, 
where theoretical gaps that allow some form of compar- 
ison can be defined. In analogy to BCS theory, one can 
then consider the state-dependent gap of the canonical 
state with the lowest quasiparticle energy (cf. Sect. Ill CP 
as a measure of pairing correlations (see, e.g., [37l. |38||): 



A, 



£ Mo = min £ M . 



(29) 



Various other prescriptions for the gap are used in the 
literature as well, in particular 



(A) = 



(30) 



which corresponds to the average of the pairing en ergy 
over the paired canonical states at the Fermi surface [25| . 
To interpret our theoretical results, we will primarily use 
the canonical gap (f2U)) . In comparison, the averaged gap 
has only slightly different values and exhibits somewhat 
smoother trends. Any exceptions to this behavior will be 
addressed explicitly in the discussion. 



E. Unitary Correlation Operator Method (UCOM) 



- (q r s(r) +s(r)q r ) 



(31) 



and 



gn = i?(r)5i 2 (r,gn) 



#(r)-((cr r q Q )(a 2 -r)+ (ovr) (a 2 -q n )) , (32) 



where 



(33) 

Qn = q - q r -_ = — { L x r - r x L ) ■ ( 34 ) 



r 1 
q r - = — (L x r r x L) 

r 2r z 



The generator g r uses the radial part of the relative mo- 
mentum operator q to create a shift in the radial direc- 
tion, while gn is constructed from the so-called orbital 
momentum, i.e., the angular part of q, and generates 
shifts perpendicular to r. The strength and range of the 
transformation is governed by the shift function s(r) and 
the tensor correlation function $(r). Rather than using 
the shift function directly, it is more practical to define 
the central correlation function R+(r) via the integral 
equation 



R+(r) ^_ 

»(0 



1. 



(35) 



which implies R+{r) ~ r + s(r) for a weakly r-dependent 
s(r). 

The unitary transformation is then written as 



C = CnC r = cxp 



i ^2 9a >i k 
j<k 



i ^ ] 9r,jk 
j<k 



(36) 

where the sum runs over all nucleon pairs. One can now 
proceed to calculate expectation values either by apply- 
ing C to the many-body state \^>) or to a given observable 
O, yielding either a correlated state or a correlated 
operator O: 

($| O I*) - (*| ClclOCnC r |$) = <*| O |$) . (37) 



The Unitary Correlation Operator Method is moti- 
vated by physical considerations on the structure and ori- 
gin of the dominant many-body correlations. The short- 
range repulsion in the central part of the NN interaction 
drives the interacting nucleon pair apart. The tensor 
interaction induces correlations between the relative dis- 
tance and the spin of the nucleon pair, leading to the 
characteristic mixing between components with relative 
orbital angular momentum L and L ± 2 in the S = 1 
channel. To imprint these correlations on a many-body 
state, we construct a unitary transformation with the 



The structure of the transformation (|36|) implies that O 
is an A-body operator in Fock space, which can be ex- 
pressed in terms of irreducible contributions for a 
specific particle number n < A via the cluster expansion 



o = c [ oc = o [1] + o [2] 



(38) 



If the range of the correlation functions is small com- 
pared to the mean inter-particle distance, we can employ 
the two-body approximation and omit negligible clus- 
ter terms beyond the second order (for details sec Rcfs. 
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@, 0, 03 )■ For the construction of the correlated Hamil- 
tonian in two-body approximation, it is then sufficient to 
consider the Hamiltonian in the two-nucleon system, 



-990 



Hint — -M r 



2/i 



V. 



(39) 



where we have already subtracted the center-of-mass ki- 
netic energy, which is not affected by the correlation pro- 
cedure. Applying the correlation operators, 



V' : 



(40) 



and collecting the two-body contributions from the cor- 
related kinetic energy and the transformed interaction, 
we obtain the effective interaction Vucom: 



Vucom = T- 



(41) 



The evaluation of the matrix elements of Vucom in a 
partial- wave basis is discussed in detail in Ref. [10( ■ 



F. SRG-generated Vucom 

In a recent pair of papers [1, [§| , we have studied the 
connection of the UCOM to the Similarity Renormaliza- 
tion Group (SRG) approach to the construction of effec- 
tive NN interactions [J]. There, the many-body Hamil- 
tonian H is evolved towards a block-diagonal structure 
in momentum space via the flow equation 



dH a 
da 



[q(a),H a ], H = H, 



where a denotes the flow parameter and 

H & = U(a)HU\a) =T int + V & . 
The dynamical generator of the flow is defined by 



77(a) = [T int , H a ] = [—,H & ] , 

2/_t 



(42) 



(43) 



(44) 



where q is the relative momentum operator. Equation 
(j4"4"]l is only appropriate for an evolution in two-body 
space, an assumption corresponding to the two-body 
approximation used in the UCOM framework. A gen- 
eralization to three-nucleon or many-nucleon forces is 
straightforward in principle, but too demanding to al- 
low the numerical evolution of realistic 3N Hamiltonians 
at present. 

In Ref. Q, we describe a procedure by which a map- 
ping between an uncorrelated trial state and the deuteron 
wave function of the SRG-evolved Hamiltonian H & de- 
fines a set of central and tensor correlation functions for 
use in the generators (|3"Tj) and (|32|) . Similar to the SRG- 
evolved V a , the resulting Vucom is uniquely determined 
by the parameter a and the parent interaction, and of- 
fers a significantly enhanced convergence. At the same 
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FIG. 1: (Color online) Ground-state energy of 120 Sn for vari- 
ous basis sizes e max (see inset) and oscillator lengths ano, us- 
ing the Gogny D1S interaction (top) and the SRG-optimized 
Vucom with a = 0.4 (bottom). 



time, however, the saturation properties are quite differ- 
ent, because Vucom does not produce the same strong 
overbinding for large nuclear masses as V & at the two- 
body level (see also Ref. [39l|). 

For practical applications, we optimize the value of at 
by considering No-Core Shell Model Calculations of the 
Tjon-line in 3 H and 4 He. For a = 0.04 fm 4 , the resulting 
ground-state energies of these nuclei are close to the ex- 
perimental values without the inclusion of 3N forces — in 
other words 3N forces generated by the UCOM transfor- 
mation more or less cancel genuine 3N forces that would 
need to supplement the parent NN interaction in order 
to obtain the correct ground-state energies. 



III. IMPLEMENTATION 
A. Basics 

Our implementation of the HFB method makes use of 
the framework established for Hartree-Fock calculations 
with Vucom 0, 0- The eigensystem (TJU) or (EU is 
solved in a spherical harmonic oscillator (SHO) configu- 
ration space, using a truncation in the oscillator quantum 
number 



e = 2n + I . 



(45) 
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where n and I denote the radial oscillator quantum 
number and orbital angular momentum, respectively. 
An original implementation of the modified Broyden's 
method discussed in Ref. [4t| is employed to accelerate 
the convergence of the HFB fields (see [4l[ for recent ap- 
plications in nuclear structure calculations). 

Figure [T] illustrates the convergence of our HFB cal- 
culations for the sample nucleus S n, using an SRG- 
generated Vucom with a = 0.4 as well as the phenomeno- 
logical Gogny D1S interaction _34j. The convergence 
rate is rather similar for both interactions as the single- 
particle basis size is increased. Since the typical energy 
gain by increasing e max from 12 to 14 (corresponding to 
13 or 15 major oscillator shells, respectively) is merely 
1—2 MeV for the Gogny D1S interaction and even smaller 
for Vucomj w e adopt the basis with e max = 12 for the re- 
mainder of this work. Using this truncation, the residual 
dependence of the ground-state energy on the oscillator 
parameter ano is already quite weak over a wide range of 
values — nevertheless, we usually carry out calculations 
for a set of oho 's to explicitly minimize the energy in this 
respect as well. 

Finally, we note that the converged ground-state en- 
ergy obtained with Vucom is in line with previous 
Hartree-Fock results [H, [Tj], providing less than half of 
the experimental binding energy. The missing binding 
energy is due to missing long-range correlations, which 
can be recovered by going beyond the mean-field ap- 
proximation Eli as well as omitted 3N or higher many- 
nucleon forces. In contrast, the Gogny D1S interaction is 
fit to experimental ground-state energies, providing the 
bulk of the 120 Sn binding energy already in a mean-field 
calculation. 



B. Intrinsic Kinetic Energy 

An interesting issue that is rarely considered in the 
literature is the effect of the center-of-mass correction 
(see, however, the detailed study in [42j|). As indicated 
in Sect. Ill At we formulate the HFB equations using the 
intrinsic kinetic energy, which can be expressed either in 
terms of two-body operators or a combination of one- and 
two-body terms: 

i<j 

The use of an intrinsic Hamiltonian has a number of con- 
sequences for our calculations. While the total ground- 
state energy is lowered by the center-of-mass correction, 
the pairing field and pairing energy obtain positive con- 
tributions from the two-body part of the Ti nt , leading 
to a reduction compared to calculations without center- 
of-mass corrections. In the plain Hartree-Fock case, the 
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FIG. 2: (Color online) Ground-state energies (top) and charge 
radii (bottom) for the tin isotopes, calculated with the Gogny 
D1S interaction: full intrinsic kinetic energy (•) and one- 
body approximation (■). Experimental values are indicated 
by black bars [M S3]. 
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FIG. 3: (Color online) Canonical single-particle spectra of tin 
isotopes: full intrinsic kinetic energy (left) and one-body ap- 
proximation (right). Calculations were done with the Gogny 
D1S interaction. 

eigenvalues of the intrinsic HF Hamiltonian can be no 
longer directly identified with single-particle energies be- 
cause its ^-dependence invalidates Koopmans' theorem. 
Indeed, if the eigenvalues are interpreted as perturbativc 
approximations to the exact separation energies, 

el F « E N+1 - E N , (48) 

one sees that the ^4-dependence of the Hamiltonian leads 
to correction terms that account for rearrangement ef- 
fects caused by the addition or removal of a nucleoli (see 
Refs. El 

In the HFB case, the need for a similar procedure be- 
comes evident in the behavior of the Lagrange parame- 
ter A, which generally assumes positive values and can 
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therefore not directly be identified with the Fermi energy 
of the system. Unfortunately, the addition or removal 
of a particle to a system with pairing is a non-trivial is- 
sue [ill [44[ , and the generalization of the aforementioned 
single-particle energy correction terms to the HFB case 
is not obvious, especially since part of the correction is 
state-dependent, whereas A is no longer associated with 
a definite single-particle level. At present, we therefore 
adopt the simpler state-independent correction discussed 
inRef. [H|, 

X co " = X - j{T int ) . (49) 

Likewise, we define corrected canonical single-particle en- 
ergies by 

z7 Y = K»-\(T int ). (50) 

In plain HF, the difference between the simple cor- 
rection and the more involved approach amounts to 
100 — 200 kcV for levels near the Fermi surface. Like- 
wise, a naive consideration of the exact separation en- 
ergy -Ejv+i — En would give rise to a correction to the 
canonical gap: 

A7 r = A^-i(T pair >, (51) 

where the kinetic pairing energy (T pa i r ) is typically 1 — 
2 MeV at most in the tin isotopes, hence this correc- 
tion would amount to 10 — 20 keV or less and has been 
omitted at present. The best strategy to avoid such con- 
ceptual uncertainties related to the perturbative defini- 
tion of the gaps and single-particle energies is the self- 
consistent calculation of odd nuclei, which will be ad- 
dressed in a subsequent publication. In addition, there 
are some general questions regarding the expectation val- 
ues of A-depcndcnt Hamiltonians in HFB states without 
sharp particle number that go beyond the scope of the 
present discussion and will be the studied elsewhere [45| . 
These caveats should be kept in mind in the following 
sections. 

To conclude our discussion, we compare our treatment 
of the intrinsic kinetic energy to the widely used one- 
body approximation, which omits the two-body contri- 
bution in Eq. l|4"Tj) altogether. Considering the tin iso- 
topic chain, the ground-state energies and charge radii 
of the two approaches differ by 1% at most (see Fig. [2]), 
while the spectroscopic structure of the resulting ground 
states is somewhat different. Figure [3] displays the canon- 
ical single-particle spectra; for the one-body approxima- 
tion, further corrections to the single-particle energies are 
typically not applied. The general trends of the single- 
particle energies are the same in both cases, but one no- 
tices that the calculations with the full intrinsic kinetic 
energy lead to a slightly reduced level density, which will 
impede pairing correlations and is expected to reduce the 
pairing energy or the gaps in comparison to the one-body 




FIG. 4: (Color online) Canonical neutron gaps of tin isotopes 
for the full intrinsic kinetic energy (•) and the one-body ap- 
proximation (W) (see text). Experimental A^ 3 '(iV) are indi- 
cated by (K) [4y|. Calculations were done with the Gogny 
D1S interaction. 

approximation. Hence, we compare the canonical neu- 
tron gaps A„ of the two approaches in Fig. [4] We find 
that the kinetic two-body term has a considerable effect, 
reducing the gap by as much as 30%, i.e. about 500 keV, 
for the mid-shell tin isotopes in calculations with the 
Gogny D1S interaction. This effect is due to the com- 
bination of the reduced level density and the repulsive 
kinetic energy contribution to the pairing field Q12p . The 
actual size of the quenching will depend on the details 
of the fit of a phcnomcnological interaction, hence one 
cannot generalize the results for Gogny D1S easily to the 
Skyrmc interactions, for instance. Comparing with the 
canonical single-particle spectrum, we also notice that 
the dips near the sub-shell closures in 106 Sn, 114 Sn, 120 Sn 
are slightly enhanced in the calculation using the full in- 
trinsic kinetic energy 

IV. Vucom AS A PAIRING FORCE 
A. Theoretical Gap Systematics 

To assess the pairing properties of Vucom, we perform 
hybrid calculations of the tin isotopic chain, using Gogny 
D1S in the particle-hole channel. Since the changes in 
ground-state energies and charge radii are minor com- 
pared to Fig. [21 we refrain from showing these results 
again, and focus directly on the canonical neutron pair- 
ing gaps. Figure O shows the pairing gaps obtained using 
Vucom for a range of parameters a = 0.03, . . . , 0.1 fm . 

The canonical gaps obtained with the various Vucom 
are about half the size of the ones obtained with Gogny 
D1S in the pairing channel (cf. Fig. 2]). In the mid-shell 
region, the experimental A^(N) are underestimated by 
about 50% as well. Varying the range of the UCOM 
transformation via <5, we find that the gaps remain very 
stable for a = 0.03 fm 4 to 0.06 fm 4 . Qualitatively, this 
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FIG. 5: (Color online) Canonical neutron gaps of tin isotopes 
for Gogny DIS+Vucom with a = 0.03 (•), 0.04 (■), 0.06 (♦), 
and 0.1 fm 4 (A). Solid lines were obtained with the full 2j n t, 
the dashed line with the one-body approximation. Experi- 
mental A( 3 \N) are indicated by (X) [isl ]. 



implies that the attractive interaction matrix elements 
that are responsible for the pairing remain mostly unaf- 
fected in this range of parameters. This is indeed the case 
for the matrix elements in the relative 1 So partial wave 
[H B El , which is expected to dominate the pairing at 
densities below saturation (48[. Moreover, the 1 Sq ma- 
trix elements hardly change at all beyond a = 0.05 fm 4 , 
and hence it is surprising that a notable reduction of 
the pairing gaps is found for a = 0.1 fm 4 , especially since 
Vucom becomes more attractive overall at the same time 
(cf. Sect. EE}. 

To understand this observation, we have to consider 
two aspects of our calculations. First, we have to 
realize that the UCOM transformation causes a pre- 
diagonalization of the two-body Hamiltonian in momen- 
tum space, focusing the attractive and repulsive strength 
of the interaction near the diagonal region that is ac- 
cessible in mean-field type calculations. Since the pair- 
ing is governed by the matrix elements near the Fermi 
surface, where the pairing tensor k is peaked, the pair- 
ing gaps are extremely sensitive to changes in the ma- 
trix elements in this very particular region of momentum 
space. The ground-state energy, however, is far less sen- 
sitive to such details. The second aspect is the formula- 
tion of the HFB method in a single-(quasi)particle basis, 
which implies that any two-nucleon state with total or- 
bital angular momentum L is a superposition of states 
with all allowed couplings of the center-of-mass and rela- 
tive orbital angular momenta. This leads to admixtures 
of relative partial waves beyond 1 Sq, which still exhibit a 
more significant a-dependence. Moreover, since the total 
spin of the nucleon pair is independent of the center-of- 
mass component of the two-body state, mostly repulsive 
spin-singlet relative partial waves are admixed to the 1 Sq 
wave, whereas the next significant contribution to pair- 
ing in nuclear matter is due to the spin-triplet 3 Pi partial 
waves (see, e.g., Ref. [H| ). 




FIG. 6: (Color online) Canonical neutron gaps of tin isotopes 
for Gogny DIS+Vucom with a = 0.04 (•), and 0.1 fm 4 (B). 

Comparison of full interaction (solid) and So partial wave 
(dashed) in the pairing channel. Experimental A' 3 '(JV) are 
indicated by (X) gf|. 

For this reason, it is instructive to consider the effect of 
a partial-wave restriction for Vucom in the pairing chan- 
nel. The canonical neutron gaps obtained using merely 
the 1 S'o matrix elements of Vucom with a — 0.04 fm 4 
and a = 0.1 fm 4 are compared to those of the full calcu- 
lations in Fig. [SI Whereas the results for the Gogny D1S 
are practically unaffected under such a restriction, we ob- 
serve a significant effect of the higher partial waves in the 
case of Vucom , which cause a decrease of the gaps by as 
much as 20%-30% in the mid-shell tin isotopes. We also 
find that the difference between the restricted and full 
calculations is enhanced for the longer-ranged correlator 
a = 0.1 fm 4 : the inclusion of the higher partial waves re- 
duces the gap by an additional 5%-10% compared to the 
calculation with a = 0.04 fm 4 . The different behavior 
of the two kinds of interactions can be explained by the 
comparably simple structure of Gogny D1S, which lacks 
q 2 , L 2 , and tensor terms that give rise to the more di- 
verse partial wave structure of realistic NN interactions. 
Moreover, in the case of Vucom, one finds that the treat- 
ment of correlations induces a host of additional tensor 
operators (see, e.g., Ref. [Toj). 

B. Comparison with SRG-evolved Interactions 

Recently, there has been an effort to use the RG- 
evolved low-momentum interaction Vjow-k as a pairing 
interaction in conjunction with the Skyrme SLy4 force 
[37l [H[ . Since this study is in the same spirit as the dis- 
cussion in this section, we have carried out similar calcu- 
lations using interactions obtained by evolving Argonnc 
VI 8 via the SRG flow equation (O to cutoffs 

A = a~ 1/4 = 1.8, 2.0, 2.4, 2.8 fm" 1 . 

While Va and Vi ow _k are slightly different conceptually, a 
Vow-k with a soft cutoff function behaves very similar to 
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FIG. 7: (Color online) Canonical neutron gaps of tin isotopes 
for Gogny DIS+Vsrg with A = 2.8 (•), 2.4 (■), 2.0 (♦), and 
1.8 fm _1 (A). Solid lines: full T int , dashed line: one-body 
approximation. Experimental A' 3 '(iV) are indicated by (X) 



an SRG-evolved interaction Moreover, their prop- 

erties with respect to binding energies in NCSM calcula- 
tions or nuclear matter approaches are similar enough to 
expect the same for their pairing properties. Compared 
to Vucom, both Viow-k and Vsrg are significantly softer, 
and require additional 3N forces to produce saturation 
in nuclear matter or finite nuclei (see Ref. [39[ and also 
Ref. i). 

The resulting canonical neutron gaps of the tin chain 
are displayed in Fig. [Jj The cutoff variation leads to 
changes on the order of 10% in the mid-shell region 
around 114 Sn. For A = 2.0 and 2.4 far 1 , the theoreti- 
cal gaps are almost identical. For these values, the flow 
affects the interaction only at momenta that are already 
decoupled from the low-energy scales relevant for nuclear 
structure. The change obtained by lowering the cutoff 
from 2.8 fm -1 can then be understood as the shift of 
repulsion to higher momenta and many-nucleon terms, 
which is typical for the SRG evolution and renders the 
interaction more attractive in the partial waves that are 
relevant for pairing. As the Vucom with a = 0.1 fin in 
the previous subsection, the corresponding SRG-evolved 
interaction with the lowest cutoff 1.8 fm -1 yields a de- 
crease of the gaps. 

The obtained results are considerably different from 
those of Ref. [38| . which presents canonical gaps A„ 
close to experimental gaps A^(N) in a series of iso- 
topic chains. Given our previous findings discussed in 
this paper, we argue that the possible reasons for this 
discrepancy are two-fold: as stated in (38|, Lesinski et 
al. use only the relative 1 Sq partial wave of Vjow-k as an 
input in the particle-particle channel at present. In the 
previous subsection, we found that higher partial waves 
can reduce the canonical gaps by as much as 20%-30% 
(although in the specific case of Vucom)- Moreover, in 
our discussion of the center-of-mass treatment in Sect. 
IIIIBI we observed a significant suppression of the gaps 



caused by the two-body term of the intrinsic kinetic en- 
ergy (cf. Fig. 2]). For this reason, we have also included 
the gaps calculated with the one-body approximation to 
the kinetic energy for A = 2.4 fm -1 in Fig. [7J These gaps 
are indeed close to the experimental values, except near 
the major shell closures where one has to include pairing 
correlations beyond the HFB approximation. The sup- 
pression of the gaps due to the kinetic two-body term 
amounts to as much as 600-700 kcV. A similar calcula- 
tion for Vucom with a = 0.04 fm 4 is included in Fig. [5] 
along with the results using the fall Ti, lt , and it exhibits 
the same effect. We note, however, that our results were 
obtained with the specific choice (|46[) for the intrinsic ki- 
netic energy, whereas Lesinski et al. used the one- plus 
two-body form (gTJ [50| • While Eqs. (gBJ) and gTJ are 
equivalent at the operator level and it has been explicitly 
shown that they lead to the same energy expectation val- 
ues in Hartree-Fock [FlJ, it is not clear that this is still 
the case in HFB calculations, in particular due to use 
of an A-dependent Hamiltonian in a state without fixed 
particle number. We will analyze this issue in detail in 
a forthcoming publication (4f|. Until then, we cannot 
rule out that three-nucleon forces or beyond mean-field 
effects like the coupling to surface vibrations may have 
important effects on the pairing gaps. 

At a first glance, our results are similar to those pre- 
sented by Barranco et al. [52|,[53[ in studies on the impact 
of particle- vibration coupling on the pairing gap in 120 Sn, 
which combine the phenomcnological Gogny D1S and 
SLy4 interactions in the particle-hole channel with the 
Argonnc V14 interaction as a pairing force. We have to 
stress, however, that this agreement might be incidental 
because there are other important aspects that need to 
be considered [HI, [HI- For instance, the use of the "bare", 
hard core AV14/18 in conjunction with phenomcnologi- 
cal forces that are essentially of low- momentum character 
in the work of Barranco et al. is certainly a consistency 
issue. 



V. FULLY SELF-CONSISTENT HFB WITH 

Vucom 

Having gained some insight on how Vucom behaves 
as a pairing force in the previous section, we now use a 
fully self-consistent HFB approach, using Vucom in the 
particle-hole as well as the particle-particle channel. 



A. Ground-State Energies and Radii 

We first consider the bulk properties of the tin iso- 
topes obtained in a fully self-consistent calculation with 
Vucom for a = 0.03, . . . , 1.0 fm 4 . Figure [5] displays the 
resulting ground-state energies and charge radii. As ex- 
pected from previous work, the nuclei are bound already 
at the mean-field level due to the explicit treatment of 
short-range correlations. The difference of 4-6 MeV per 
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FIG. 8: (Color online) Ground-state energies (top) and charge 
radii (bottom) for the tin isotopes, calculated for Vucom 
with a = 0.03 (•), 0.04 (■), 0.05 (♦), 0.06 (A), 0.08 (OV and 
0.1 fm 4 (□). Experimental values ( ) taken from [4q.l47j. 



nucleon from experimental data is due to long-range cor- 
relations that are not described by the UCOM correla- 
tion operators and can be described by beyond mean-field 
methods like many-body perturbation theory, as demon- 
strated successfully in [IE, E3 . The increase of the bind- 
ing energy with a implies that longer-ranged correlations 
are shifted from the many-body state into the correlation 
operators; it roughly corresponds to the cutoff depen- 
dence of results obtained with pure two-body Viow-k or 
SRG-evolved interactions. 

In contrast to Ref. [IE], however, in which the UCOM 
transformation was constructed by a different strategy, 
the new SRG-generated correlation functions provide a 
significant improvement of the charge radii, which lie 
within roughly 10% of experimental data and correctly 
reproduce the experimentally observed trends over the 
isotopic chain (see also @). Since the radius is a long- 
ranged operator, it is rather insensitive to variations of 
a, as evident from Fig. [5] (note the scale of the plot). 



B. Gaps 

In Fig. [9l we show the canonical neutron gaps of 
the tin isotopes. As in the hybrid calculations in Sect. 
IV! the gaps are very stable under variations of a. It 
is noteworthy that this is even the case for the long- 
ranged Vucom with a = 0.1 fm 4 . This improved stability 
could be a signal of the improved consistency, because 
unlike in the hybrid Gogny DIS+Vucom calculations, 
the single-particle spectrum is directly affected by the 
variation of a as well. Compared to the experimentally 
extracted A^ 3 '(iV), the theoretical gaps are significantly 
lower, ranging from below 100 keV from the outer tin 
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FIG. 9: (Color online) Canonical gaps in 

the tin isotopes for Vucom with a = 

0.03 (•), 0.04 (■), 0.05 (♦), 0.06 (a), 0.08 (O), and 

0.1 fm 4 (□), compared to experimental A (3) (iV) (*) 



isotopes to roughly 400 keV in the mid-shell region. In 
addition, we find clear signals of each sub-shell closure in 
the tin nuclei. 

These findings can be understood if we consider the 
canonical single-particle spectra, which are shown for 
Vucom with a = 0.04 fm 4 in Fig. [I0l We note that the 
canonical neutron energies are spread over an interval on 
the order of 10 MeV, which is about twice as large as for 
the purely phcnomcnological calculations with the Gogny 
D1S interaction in Fig. [3] The discrepancy between our 
spectra and the experimentally extracted single-particle 
levels of 132 Sn, which are included in Fig. [10] for refer- 
ence, is even more severe. Such a low level density is 
a common feature of soft NN interactions, and consis- 
tent with previous studies using Vucom EE EE E3 ■ 
From BCS theory, it is well-known that the formation of 
Cooper pairs strongly depends on a sufficiently high level 
density in the region of the Fermi surface. For the BCS 
gap, this is reflected by the relation [56[ 



A 



1 



cxp - 



IsKo) 



(52) 



where g is the strength of the (attractive) pairing interac- 
tion and n(0) the level density at the Fermi surface. Con- 
sequently, we see that the low level density generated by 
Vucom presents a major obstacle to nuclear pairing. To 
obtain more realistic single-particle spectra, we will have 
to account for long-range correlations that are presently 
not described by either the correlation operators or the 
relatively simple many-body space, as well as three- or 
possibly higher many-nucleon forces. 

Comparing the theoretical and experimental single- 
particle levels for 132 Sn, we see that the 0g 7 / 2 and 0hn/ 2 
shells show a particularly large deviation, whereas the re- 
maining levels and their splittings are reproduced rather 
well. Comparing with Fig. [5J we see that for these levels 
the canonical gaps are strongly suppressed as well, com- 
pared to the mid-shell region where the lowest canonical 
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FIG. 10: (Color online) Corrected canonical neutron single- 
particle energies in the tin isotopes, calculated for a Vucom 
with a = 0.04 fm 4 . Experimental single-particle energies of 
132 Sn ("Exp.") are included for reference [53]. 



quasiparticlc energies are associated with s or d orbitals. 
This strong dependence on the single-particle angular 
momenta suggests a significant influence of the tensor 
structure of Vucom- 

To gain further insight, we first compare the canonical 
neutron gaps to the average gaps defined by Eq. ([50)) in 
Fig. [TTJ We find that the latter are practically constant 
over the tin isotopic chain, suggesting that the shells with 
low single-particle j provide the essential contribution to 
this quantity, except at the sub-shell closures, where the 
pairing collapses and the solution is reduced to the HF 
case. Next, we revisit the restriction of the pairing in- 
teraction to the relative 1 Sq partial wave, which was dis- 
cussed for the hybrid calculation in Sect. IIVI This elim- 
inates tensor effects in the pairing field, and as a result, 
we find a substantial increase in the canonical gaps in 
Fig. QT] Considering that the same restriction of Vucom 
in the Gogny DIS+Vjcom calculation presented in Fig. 
[5] only caused a much smaller increase of the gap, we 
have to conclude that the interplay with the ph interac- 
tion via self-consistency effects plays an important role as 
well. Furthermore, we see that without the tensor inter- 
action, the canonical and average gaps arc very similar, 
suggesting that the pairing is balanced more uniformly 
over all shells. 

Finally, we have tested the sensitivity of these results 
to changes in a. In Fig. [TTJ we show only the calculation 
using the long-ranged Vucom with a = 0.1 fm 4 for com- 
parison. As expected from the previous discussions, we 
find only minor differences due to the a-variation. The 
most notable changes occur in the high-j subshells, and 
are in line with our previous statement that higher-lying 
partial waves require larger a to become stable under a 
variations (if at all) . The negligible difference in both the 
canonical and average gaps when the pairing interaction 
is restricted to the relative 1 Sq waves serves as a further 
confirmation. 
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FIG. 11: (Color online) Canonical (top) and average gaps 
(bottom) of tin isotopes for Vucom with a = 0.04 (•), and 
0.1 fm 4 (H). Comparison of full interaction (solid) and 1 So 
partial wave (dashed). Experimental A' 3 '(iV) are indicated 
by(*) 0. 



The discussed observations underline the significance 
of the tensor structure of a realistic (albeit effective) NN 
interaction for the pairing correlations in finite nuclei. Al- 
though the HFB approach considers these effects merely 
on a mean-field level, and the possibility of a differ- 
ent behavior in more refined many-body methods exists, 
it should nevertheless be clear that the rather simplis- 
tic phenomenological forces obscure aspects of the two- 
nucleon physics that may prove to be very important for 
spectroscopic observables. 



VI. CONCLUSIONS 

In this paper, we have presented a fully self-consistent 
HFB approach based on an intrinsic Hamiltonian. By 
using effective NN interactions derived from the realis- 
tic Argonnc V18 interaction, we are able to maintain a 
stringent link to low-energy observables of the strong in- 
teraction, i.e., NN scattering phase shifts and deuteron 
properties. Since the same interactions can also be used 
in other Hartree-Fock based approaches and more re- 
fined ab-initio many-body methods like the No Core Shell 
Model [Hi or the Coupled Cluster Method [H, this 
opens important perspectives for cross-checking nuclear 
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structure studies. 

Focusing on the tin isotopic chain, we have calculated 
theoretical pairing gaps, and studied the effects of com- 
monly used approximations in the center-of-mass treat- 
ment as well as the ansatz for phenomenological pairing 
interactions in Density Functional Theory [25[ . We have 
discussed aspects of the tensorial structure of realistic NN 
interactions whose impact on the spectroscopic proper- 
ties of finite nuclei is expected to be significant and that 
are described only inadequately by existing phenomeno- 
logical functionals of the Skyrme or Gogny type [H, [13] . 
In particular, we have demonstrated that the inclusion 
of all partial waves of the NN interaction as well as the 
repulsive contribution of the intrinsic kinetic energy in 
the pairing field have a significant effect on the gaps. In 
the latter case, however, it remains to be seen whether 
this is affected by the particular choice of the intrinsic 
kinetic energy operator (50l |. and on the treatment of its 
A dependence [451 ] . 

In a fully self-consistent approach using Vucom in the 
particle-hole and the particle-particle channel, the low 
density of single-particle levels near the Fermi surface 
proves to be a strong impediment to nucleoli Cooper 
pairing. This low level density, implying a low effective 
mass as well, is a general feature of soft, non-local in- 
teractions, and consistent with previous studies (l3l. [l5|. 
Consequently, we conclude that beyond mean-field effects 
like the coupling to surface vibrations indeed play an im- 
portant role in nuclear pairing. Such a coupling would 
lead to a dressing of the single-particle energies, and is 
expected to improve the level density near the Fermi sur- 
face, which would at least partially overcome the effects 
of the non- locality of Vucom and similar interactions. 

In principle, there are two directions in which an ex- 
tension of our framework can proceed: the inclusion of 
higher many-body forces, and the use of a more sophis- 
ticated many-body method. One of the ultimate aims 
of effective interaction methods is to obtain results that 
are independent of the control parameters of the trans- 
formation, in our case a. In this sense, the transforma- 
tion merely yields a unitarily equivalent representation of 
low-energy QCD that is more suitable to the used many- 
body method. To achieve true consistency, however, one 
would need (i) to start from a consistent set of NN and 
higher nucleonic interactions, and (ii) to include them 
in the unitary transformation of the many-body states 
or operators. The former requirement can be met by 
using N3LO interactions derived in chiral effective field 
theory. While the required consistent set of NN, 3N, 
and 4N interactions have been worked out in principle 
(see, e.g., [5!l[6(|), the complex structure of the full 3N 
terms and the computational demands for handling a 4N 
interaction have thus far prevented their use in many- 
body calculations. The second aspect complicates mat- 
ters even further, because the inclusion of the 3N force 
in either UCOM or SRG transformations is a formidable 
challenge. At present, one therefore hopes that such a 
transformation may render the 3N interaction less impor- 



tant in actual calculations, so that a simpler model like 
the chiral N2LO interaction may reproduce the required 
effects after a readjustment of its parameters [6l[. In this 
spirit, we have studied the use of a regularized 3N con- 
tact force in conjunction with Vucom as a first step [l4[. 
While initial HFB results including 3N forces together 
with first-generation UCOM interactions are available in 
[62| . we are preparing a paper on such calculations with 
the new SRC-generated correlation functions discussed 

in sect, inn 

As for the many-body methods, the use of more refined 
approximations improves the Hilbert space, enabling it 
to describe residual long-range correlations that are not 
explicitly treated by the UCOM correlation operators. 
In the context of HFB, a straightforward extension is 
the use of projection techniques to restore symmetries 
that are spontaneously broken in the calculated ground 
state. The simplest example is Particle Number Projec- 
tion, which can be implemented rather easily, because 
the general structure of the HFB eigenvalue problem is 
preserved (see, e.g., [63j). Another approach that can de- 
scribe additional correlations as well as collective behav- 
ior is the Quasiparticlc Random Phase Approximation, 
using the HFB ground state as a starting point. Both 
of these methods were explored for the first-generation 
UCOM interactions in [631, an d wm be the subject of 
studies using the SRG-generated Vucom interactions in 
subsequent papers. 
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APPENDIX A: EXPRESSIONS FOR THE 
GOGNY INTERACTION 

We briefly recall the parameterization of the Gogny 
interactions, which is given by [64[ 

V12 = ^exp {W l + BiP a - H t P T - MiP^Pr) 

i=i V / 

+ t (l + x Pa)p a (R) S 3 (r) 

+ iW LS ((Ti + (T2)- VX<y< 3 >(r) V, (Al) 

where P a and P T are the usual spin and isospin exchange 
operators, e.g., 

P a = \{l + <T 1 -a 2 ), (A2) 
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and the relative and center-of-mass positions are given 
by 

r = n-r 2 , (A3) 
R=i(r 1 +r 2 ), (A4) 

and the gradient operator in the relative coordinates is 

V = Vi - V 2 . (A5) 



with 



I, 



Til 2iTl2^2j^i2i n 2^2 



drr 2 p a (r)R nill (r)R no j 2 (r)R nlilli (r)R nl2V2 (r^ . (A9) 



The angular-momentum and isospin projections have 
been suppressed since the matrix element does not de- 
pend on them. 



1. Two-Body Matrix Elements 



2. Fields 



In Ref. [1 31 ] . we have provided expressions for 
the evaluation of two-body matrix elements in a rel- 
ative LS'-couplcd spherical harmonic oscillator basis 
\n(LS)J MTMt) ■, and the subsequent Talmi transfor- 
mation to obtain jj-coupled matrix elements for use in 
Hartree-Fock, Hartree-Fock-Bogoliubov and their exten- 
sions. With these expressions, the evaluation of two-body 
matrix elements of the finite range terms in Eq. (|A1[) is 
straightforward. 

While the spin-orbit interaction in (IA1[) is given in a 
form which facilitates the calculation of particle-hole and 
particle-particle fields for use in self-consistent field meth- 
ods, the calculation of the corresponding relative two- 
body matrix element is somewhat more involved than for 
the usual relative (L ■ S) interaction. Suppressing isospin 
indices as well as the angular-momentum projection M, 
and using rotational invariance, we obtain 



iW LS (n(LS)J\ [ V x 
9 



S(r) 



V • (o-i + <x 2 ) \n'(L'S') J) 



x N nl N nll C 3 J 2 (0)£. 3 n / , 2 (0)5 L1 5 LLI S sl , (A6) 

where /^^(cc) are the Laguerre polynomials, auo is the 
oscillator length of the relative basis, and 



2n\ 



a! o r(n + L + f) 



(A7) 



The density-dependent matrix clement is most conve- 
niently evaluated in a jj-coupled basis, and one finds 

(n 1 l 1 j 1 ,n 2 l2j2;JT\ v[p] \n' 1 l' 1 j[,n' 2 l' 2 j' 2 ; JT) 
1 T M^hKj'iK 

— 2^ + x °> 4 7r (2J+ 1) ^ih^h-y^n'^ 



x (7i5i2-5l J0 >0'i5J2-3l J0 > 

+ (l + (-If) (ji\h\\Jl){j[U\\Jl)\ , (A8) 



For the sake of efficiency, we calculate the particlc-holc 
and particle-particle fields of the density-dependent in- 
teraction as in density functional approaches rather than 
by contracting p or k with the matrix element (|A8[) (see, 
eg-, Ref. 0). Using 



(A10) 



kk' 



K r( r ) ^Y. K Vk'^r{r)^- klT {r) 1 (All) 



kk' 



where k denotes a time-reversed state and r = p,n, the 
matrix elements of the fields are given by 



(A12) 



and 



Here. 



Tl k , = / d d r^ T (r-)r T (r)^ T (r-) 



A T kk , = / d 3 r A T (r)^ kT (r)^ k ' T (r) . 



(A13) 



f T (r) =t 



1 + yj P a+1 (r) - \x + -j p a (r)p T (r) 
(1 - x ) p<*-\r) L\r) + ]T K<(r)| 2 ) 



+ t a \x + - j p a - L (r)p p (r)p n (r) , (A14) 

where the parts proportional to a constitute the rear- 
rangement term due to the density dependence of the 
interaction, and 

A T (r) = h (l-x )p a (r) KT (r). (A15) 

In the case of spherical symmetry, the densities are 
reduced to 



Ijnn' 



rljnn' 



-^-Pnn' R nl{r)Rn'l{r) , 

gj + l 

47T 



(A16) 



{-l) l K l Z,R nl {r)R n>l {r), (A17) 
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where the phase (—1)' appears due to using the proper- 
ties of the spherical harmonics under time reversal, and 
Rni (r) are radial spherical harmonic oscillator wavefunc- 
tions. Likewise, the fields are 



r^ = J drR nl {r)f T {r)R n n{r), 
= J drR nl {r)KP{r)R nll {r) 

with 

f T (r) =t 



(A18) 
(A19) 



(l + y) P a+ \r) - (x + iy a (r) P r(r) 
t ja(l-x )p a -\r) (W)+|>t'M) 



t a [ x + | ) p a 1 {r)pp{r)p n {r) 



(A20) 



where we have used that p(r) and n(r) are real, and 



AW)(r) = i*o (1 - so) p Q (r)(-l)^ r (r) . (A21) 
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